3.2.30 \(\int \frac {\sec ^4(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx\) [130]

3.2.30.1 Optimal result
3.2.30.2 Mathematica [C] (verified)
3.2.30.3 Rubi [A] (verified)
3.2.30.4 Maple [B] (verified)
3.2.30.5 Fricas [A] (verification not implemented)
3.2.30.6 Sympy [F]
3.2.30.7 Maxima [F(-1)]
3.2.30.8 Giac [F(-2)]
3.2.30.9 Mupad [F(-1)]

3.2.30.1 Optimal result

Integrand size = 23, antiderivative size = 181 \[ \int \frac {\sec ^4(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=-\frac {9 \text {arctanh}\left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {a+a \cos (c+d x)}}\right )}{8 \sqrt {a} d}+\frac {\sqrt {2} \text {arctanh}\left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {2} \sqrt {a+a \cos (c+d x)}}\right )}{\sqrt {a} d}+\frac {7 \tan (c+d x)}{8 d \sqrt {a+a \cos (c+d x)}}-\frac {\sec (c+d x) \tan (c+d x)}{12 d \sqrt {a+a \cos (c+d x)}}+\frac {\sec ^2(c+d x) \tan (c+d x)}{3 d \sqrt {a+a \cos (c+d x)}} \]

output
-9/8*arctanh(sin(d*x+c)*a^(1/2)/(a+a*cos(d*x+c))^(1/2))/d/a^(1/2)+arctanh( 
1/2*sin(d*x+c)*a^(1/2)*2^(1/2)/(a+a*cos(d*x+c))^(1/2))*2^(1/2)/d/a^(1/2)+7 
/8*tan(d*x+c)/d/(a+a*cos(d*x+c))^(1/2)-1/12*sec(d*x+c)*tan(d*x+c)/d/(a+a*c 
os(d*x+c))^(1/2)+1/3*sec(d*x+c)^2*tan(d*x+c)/d/(a+a*cos(d*x+c))^(1/2)
 
3.2.30.2 Mathematica [C] (verified)

Result contains complex when optimal does not.

Time = 1.42 (sec) , antiderivative size = 385, normalized size of antiderivative = 2.13 \[ \int \frac {\sec ^4(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\frac {\cos \left (\frac {1}{2} (c+d x)\right ) \sec ^3(c+d x) \left (9 \cos (c+d x) \left (9 \sqrt {2} \log \left (i-\sqrt {2} e^{\frac {1}{2} i (c+d x)}-i e^{i (c+d x)}\right )-9 \sqrt {2} \log \left (i+\sqrt {2} e^{\frac {1}{2} i (c+d x)}-i e^{i (c+d x)}\right )-32 \log \left (\cos \left (\frac {1}{4} (c+d x)\right )-\sin \left (\frac {1}{4} (c+d x)\right )\right )+32 \log \left (\cos \left (\frac {1}{4} (c+d x)\right )+\sin \left (\frac {1}{4} (c+d x)\right )\right )\right )+3 \cos (3 (c+d x)) \left (9 \sqrt {2} \log \left (i-\sqrt {2} e^{\frac {1}{2} i (c+d x)}-i e^{i (c+d x)}\right )-9 \sqrt {2} \log \left (i+\sqrt {2} e^{\frac {1}{2} i (c+d x)}-i e^{i (c+d x)}\right )-32 \log \left (\cos \left (\frac {1}{4} (c+d x)\right )-\sin \left (\frac {1}{4} (c+d x)\right )\right )+32 \log \left (\cos \left (\frac {1}{4} (c+d x)\right )+\sin \left (\frac {1}{4} (c+d x)\right )\right )\right )+4 \left (78 \sin \left (\frac {1}{2} (c+d x)\right )-25 \sin \left (\frac {3}{2} (c+d x)\right )+21 \sin \left (\frac {5}{2} (c+d x)\right )\right )\right )}{192 d \sqrt {a (1+\cos (c+d x))}} \]

input
Integrate[Sec[c + d*x]^4/Sqrt[a + a*Cos[c + d*x]],x]
 
output
(Cos[(c + d*x)/2]*Sec[c + d*x]^3*(9*Cos[c + d*x]*(9*Sqrt[2]*Log[I - Sqrt[2 
]*E^((I/2)*(c + d*x)) - I*E^(I*(c + d*x))] - 9*Sqrt[2]*Log[I + Sqrt[2]*E^( 
(I/2)*(c + d*x)) - I*E^(I*(c + d*x))] - 32*Log[Cos[(c + d*x)/4] - Sin[(c + 
 d*x)/4]] + 32*Log[Cos[(c + d*x)/4] + Sin[(c + d*x)/4]]) + 3*Cos[3*(c + d* 
x)]*(9*Sqrt[2]*Log[I - Sqrt[2]*E^((I/2)*(c + d*x)) - I*E^(I*(c + d*x))] - 
9*Sqrt[2]*Log[I + Sqrt[2]*E^((I/2)*(c + d*x)) - I*E^(I*(c + d*x))] - 32*Lo 
g[Cos[(c + d*x)/4] - Sin[(c + d*x)/4]] + 32*Log[Cos[(c + d*x)/4] + Sin[(c 
+ d*x)/4]]) + 4*(78*Sin[(c + d*x)/2] - 25*Sin[(3*(c + d*x))/2] + 21*Sin[(5 
*(c + d*x))/2])))/(192*d*Sqrt[a*(1 + Cos[c + d*x])])
 
3.2.30.3 Rubi [A] (verified)

Time = 1.18 (sec) , antiderivative size = 206, normalized size of antiderivative = 1.14, number of steps used = 16, number of rules used = 15, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.652, Rules used = {3042, 3258, 3042, 3463, 27, 3042, 3463, 27, 3042, 3464, 3042, 3128, 219, 3252, 219}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\sec ^4(c+d x)}{\sqrt {a \cos (c+d x)+a}} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {1}{\sin \left (c+d x+\frac {\pi }{2}\right )^4 \sqrt {a \sin \left (c+d x+\frac {\pi }{2}\right )+a}}dx\)

\(\Big \downarrow \) 3258

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\int \frac {(a-5 a \cos (c+d x)) \sec ^3(c+d x)}{\sqrt {\cos (c+d x) a+a}}dx}{6 a}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\int \frac {a-5 a \sin \left (c+d x+\frac {\pi }{2}\right )}{\sin \left (c+d x+\frac {\pi }{2}\right )^3 \sqrt {\sin \left (c+d x+\frac {\pi }{2}\right ) a+a}}dx}{6 a}\)

\(\Big \downarrow \) 3463

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {\int -\frac {3 \left (7 a^2-a^2 \cos (c+d x)\right ) \sec ^2(c+d x)}{2 \sqrt {\cos (c+d x) a+a}}dx}{2 a}+\frac {a \tan (c+d x) \sec (c+d x)}{2 d \sqrt {a \cos (c+d x)+a}}}{6 a}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {a \tan (c+d x) \sec (c+d x)}{2 d \sqrt {a \cos (c+d x)+a}}-\frac {3 \int \frac {\left (7 a^2-a^2 \cos (c+d x)\right ) \sec ^2(c+d x)}{\sqrt {\cos (c+d x) a+a}}dx}{4 a}}{6 a}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {a \tan (c+d x) \sec (c+d x)}{2 d \sqrt {a \cos (c+d x)+a}}-\frac {3 \int \frac {7 a^2-a^2 \sin \left (c+d x+\frac {\pi }{2}\right )}{\sin \left (c+d x+\frac {\pi }{2}\right )^2 \sqrt {\sin \left (c+d x+\frac {\pi }{2}\right ) a+a}}dx}{4 a}}{6 a}\)

\(\Big \downarrow \) 3463

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {a \tan (c+d x) \sec (c+d x)}{2 d \sqrt {a \cos (c+d x)+a}}-\frac {3 \left (\frac {\int -\frac {\left (9 a^3-7 a^3 \cos (c+d x)\right ) \sec (c+d x)}{2 \sqrt {\cos (c+d x) a+a}}dx}{a}+\frac {7 a^2 \tan (c+d x)}{d \sqrt {a \cos (c+d x)+a}}\right )}{4 a}}{6 a}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {a \tan (c+d x) \sec (c+d x)}{2 d \sqrt {a \cos (c+d x)+a}}-\frac {3 \left (\frac {7 a^2 \tan (c+d x)}{d \sqrt {a \cos (c+d x)+a}}-\frac {\int \frac {\left (9 a^3-7 a^3 \cos (c+d x)\right ) \sec (c+d x)}{\sqrt {\cos (c+d x) a+a}}dx}{2 a}\right )}{4 a}}{6 a}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {a \tan (c+d x) \sec (c+d x)}{2 d \sqrt {a \cos (c+d x)+a}}-\frac {3 \left (\frac {7 a^2 \tan (c+d x)}{d \sqrt {a \cos (c+d x)+a}}-\frac {\int \frac {9 a^3-7 a^3 \sin \left (c+d x+\frac {\pi }{2}\right )}{\sin \left (c+d x+\frac {\pi }{2}\right ) \sqrt {\sin \left (c+d x+\frac {\pi }{2}\right ) a+a}}dx}{2 a}\right )}{4 a}}{6 a}\)

\(\Big \downarrow \) 3464

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {a \tan (c+d x) \sec (c+d x)}{2 d \sqrt {a \cos (c+d x)+a}}-\frac {3 \left (\frac {7 a^2 \tan (c+d x)}{d \sqrt {a \cos (c+d x)+a}}-\frac {9 a^2 \int \sqrt {\cos (c+d x) a+a} \sec (c+d x)dx-16 a^3 \int \frac {1}{\sqrt {\cos (c+d x) a+a}}dx}{2 a}\right )}{4 a}}{6 a}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {a \tan (c+d x) \sec (c+d x)}{2 d \sqrt {a \cos (c+d x)+a}}-\frac {3 \left (\frac {7 a^2 \tan (c+d x)}{d \sqrt {a \cos (c+d x)+a}}-\frac {9 a^2 \int \frac {\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right ) a+a}}{\sin \left (c+d x+\frac {\pi }{2}\right )}dx-16 a^3 \int \frac {1}{\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right ) a+a}}dx}{2 a}\right )}{4 a}}{6 a}\)

\(\Big \downarrow \) 3128

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {a \tan (c+d x) \sec (c+d x)}{2 d \sqrt {a \cos (c+d x)+a}}-\frac {3 \left (\frac {7 a^2 \tan (c+d x)}{d \sqrt {a \cos (c+d x)+a}}-\frac {9 a^2 \int \frac {\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right ) a+a}}{\sin \left (c+d x+\frac {\pi }{2}\right )}dx+\frac {32 a^3 \int \frac {1}{2 a-\frac {a^2 \sin ^2(c+d x)}{\cos (c+d x) a+a}}d\left (-\frac {a \sin (c+d x)}{\sqrt {\cos (c+d x) a+a}}\right )}{d}}{2 a}\right )}{4 a}}{6 a}\)

\(\Big \downarrow \) 219

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {a \tan (c+d x) \sec (c+d x)}{2 d \sqrt {a \cos (c+d x)+a}}-\frac {3 \left (\frac {7 a^2 \tan (c+d x)}{d \sqrt {a \cos (c+d x)+a}}-\frac {9 a^2 \int \frac {\sqrt {\sin \left (c+d x+\frac {\pi }{2}\right ) a+a}}{\sin \left (c+d x+\frac {\pi }{2}\right )}dx-\frac {16 \sqrt {2} a^{5/2} \text {arctanh}\left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {2} \sqrt {a \cos (c+d x)+a}}\right )}{d}}{2 a}\right )}{4 a}}{6 a}\)

\(\Big \downarrow \) 3252

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {a \tan (c+d x) \sec (c+d x)}{2 d \sqrt {a \cos (c+d x)+a}}-\frac {3 \left (\frac {7 a^2 \tan (c+d x)}{d \sqrt {a \cos (c+d x)+a}}-\frac {-\frac {18 a^3 \int \frac {1}{a-\frac {a^2 \sin ^2(c+d x)}{\cos (c+d x) a+a}}d\left (-\frac {a \sin (c+d x)}{\sqrt {\cos (c+d x) a+a}}\right )}{d}-\frac {16 \sqrt {2} a^{5/2} \text {arctanh}\left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {2} \sqrt {a \cos (c+d x)+a}}\right )}{d}}{2 a}\right )}{4 a}}{6 a}\)

\(\Big \downarrow \) 219

\(\displaystyle \frac {\tan (c+d x) \sec ^2(c+d x)}{3 d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {a \tan (c+d x) \sec (c+d x)}{2 d \sqrt {a \cos (c+d x)+a}}-\frac {3 \left (\frac {7 a^2 \tan (c+d x)}{d \sqrt {a \cos (c+d x)+a}}-\frac {\frac {18 a^{5/2} \text {arctanh}\left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {a \cos (c+d x)+a}}\right )}{d}-\frac {16 \sqrt {2} a^{5/2} \text {arctanh}\left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {2} \sqrt {a \cos (c+d x)+a}}\right )}{d}}{2 a}\right )}{4 a}}{6 a}\)

input
Int[Sec[c + d*x]^4/Sqrt[a + a*Cos[c + d*x]],x]
 
output
(Sec[c + d*x]^2*Tan[c + d*x])/(3*d*Sqrt[a + a*Cos[c + d*x]]) - ((a*Sec[c + 
 d*x]*Tan[c + d*x])/(2*d*Sqrt[a + a*Cos[c + d*x]]) - (3*(-1/2*((18*a^(5/2) 
*ArcTanh[(Sqrt[a]*Sin[c + d*x])/Sqrt[a + a*Cos[c + d*x]]])/d - (16*Sqrt[2] 
*a^(5/2)*ArcTanh[(Sqrt[a]*Sin[c + d*x])/(Sqrt[2]*Sqrt[a + a*Cos[c + d*x]]) 
])/d)/a + (7*a^2*Tan[c + d*x])/(d*Sqrt[a + a*Cos[c + d*x]])))/(4*a))/(6*a)
 

3.2.30.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 219
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))* 
ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x] /; FreeQ[{a, b}, x] && NegQ[a/b] && (Gt 
Q[a, 0] || LtQ[b, 0])
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3128
Int[1/Sqrt[(a_) + (b_.)*sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[-2/d 
Subst[Int[1/(2*a - x^2), x], x, b*(Cos[c + d*x]/Sqrt[a + b*Sin[c + d*x]])], 
 x] /; FreeQ[{a, b, c, d}, x] && EqQ[a^2 - b^2, 0]
 

rule 3252
Int[Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]/((c_.) + (d_.)*sin[(e_.) + ( 
f_.)*(x_)]), x_Symbol] :> Simp[-2*(b/f)   Subst[Int[1/(b*c + a*d - d*x^2), 
x], x, b*(Cos[e + f*x]/Sqrt[a + b*Sin[e + f*x]])], x] /; FreeQ[{a, b, c, d, 
 e, f}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0]
 

rule 3258
Int[((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_)/Sqrt[(a_) + (b_.)*sin[(e_. 
) + (f_.)*(x_)]], x_Symbol] :> Simp[(-d)*Cos[e + f*x]*((c + d*Sin[e + f*x]) 
^(n + 1)/(f*(n + 1)*(c^2 - d^2)*Sqrt[a + b*Sin[e + f*x]])), x] - Simp[1/(2* 
b*(n + 1)*(c^2 - d^2))   Int[(c + d*Sin[e + f*x])^(n + 1)*(Simp[a*d - 2*b*c 
*(n + 1) + b*d*(2*n + 3)*Sin[e + f*x], x]/Sqrt[a + b*Sin[e + f*x]]), x], x] 
 /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] 
&& NeQ[c^2 - d^2, 0] && LtQ[n, -1] && IntegerQ[2*n]
 

rule 3463
Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((A_.) + (B_.)*sin[(e_.) + 
(f_.)*(x_)])*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Sim 
p[(B*c - A*d)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x])^(n 
+ 1)/(f*(n + 1)*(c^2 - d^2))), x] + Simp[1/(b*(n + 1)*(c^2 - d^2))   Int[(a 
 + b*Sin[e + f*x])^m*(c + d*Sin[e + f*x])^(n + 1)*Simp[A*(a*d*m + b*c*(n + 
1)) - B*(a*c*m + b*d*(n + 1)) + b*(B*c - A*d)*(m + n + 2)*Sin[e + f*x], x], 
 x], x] /; FreeQ[{a, b, c, d, e, f, A, B, m}, x] && NeQ[b*c - a*d, 0] && Eq 
Q[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0] && LtQ[n, -1] && (IntegerQ[n] || EqQ[m 
 + 1/2, 0])
 

rule 3464
Int[((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])/(Sqrt[(a_) + (b_.)*sin[(e_.) + 
(f_.)*(x_)]]*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])), x_Symbol] :> Simp[(A 
*b - a*B)/(b*c - a*d)   Int[1/Sqrt[a + b*Sin[e + f*x]], x], x] + Simp[(B*c 
- A*d)/(b*c - a*d)   Int[Sqrt[a + b*Sin[e + f*x]]/(c + d*Sin[e + f*x]), x], 
 x] /; FreeQ[{a, b, c, d, e, f, A, B}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - 
 b^2, 0] && NeQ[c^2 - d^2, 0]
 
3.2.30.4 Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(882\) vs. \(2(152)=304\).

Time = 2.04 (sec) , antiderivative size = 883, normalized size of antiderivative = 4.88

method result size
default \(\text {Expression too large to display}\) \(883\)

input
int(sec(d*x+c)^4/(a+cos(d*x+c)*a)^(1/2),x,method=_RETURNVERBOSE)
 
output
1/6*cos(1/2*d*x+1/2*c)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)*(-24*a*(16*2^(1/2)*l 
n(4*(a^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)+a)/cos(1/2*d*x+1/2*c))-9*ln(4/ 
(2*cos(1/2*d*x+1/2*c)+2^(1/2))*(2^(1/2)*a*cos(1/2*d*x+1/2*c)+2^(1/2)*(a*si 
n(1/2*d*x+1/2*c)^2)^(1/2)*a^(1/2)+2*a))-9*ln(-4/(2*cos(1/2*d*x+1/2*c)-2^(1 
/2))*(2^(1/2)*a*cos(1/2*d*x+1/2*c)-2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)* 
a^(1/2)-2*a)))*sin(1/2*d*x+1/2*c)^6+(576*2^(1/2)*ln(4*(a^(1/2)*(a*sin(1/2* 
d*x+1/2*c)^2)^(1/2)+a)/cos(1/2*d*x+1/2*c))*a-324*ln(4/(2*cos(1/2*d*x+1/2*c 
)+2^(1/2))*(2^(1/2)*a*cos(1/2*d*x+1/2*c)+2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^ 
(1/2)*a^(1/2)+2*a))*a-324*ln(-4/(2*cos(1/2*d*x+1/2*c)-2^(1/2))*(2^(1/2)*a* 
cos(1/2*d*x+1/2*c)-2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)*a^(1/2)-2*a))*a+ 
168*2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)*a^(1/2))*sin(1/2*d*x+1/2*c)^4+( 
-288*2^(1/2)*ln(4*(a^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)+a)/cos(1/2*d*x+1 
/2*c))*a+162*ln(4/(2*cos(1/2*d*x+1/2*c)+2^(1/2))*(2^(1/2)*a*cos(1/2*d*x+1/ 
2*c)+2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1/2)*a^(1/2)+2*a))*a+162*ln(-4/(2*c 
os(1/2*d*x+1/2*c)-2^(1/2))*(2^(1/2)*a*cos(1/2*d*x+1/2*c)-2^(1/2)*(a*sin(1/ 
2*d*x+1/2*c)^2)^(1/2)*a^(1/2)-2*a))*a-160*2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2) 
^(1/2)*a^(1/2))*sin(1/2*d*x+1/2*c)^2+48*2^(1/2)*ln(4*(a^(1/2)*(a*sin(1/2*d 
*x+1/2*c)^2)^(1/2)+a)/cos(1/2*d*x+1/2*c))*a-27*ln(4/(2*cos(1/2*d*x+1/2*c)+ 
2^(1/2))*(2^(1/2)*a*cos(1/2*d*x+1/2*c)+2^(1/2)*(a*sin(1/2*d*x+1/2*c)^2)^(1 
/2)*a^(1/2)+2*a))*a-27*ln(-4/(2*cos(1/2*d*x+1/2*c)-2^(1/2))*(2^(1/2)*a*...
 
3.2.30.5 Fricas [A] (verification not implemented)

Time = 0.28 (sec) , antiderivative size = 263, normalized size of antiderivative = 1.45 \[ \int \frac {\sec ^4(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\frac {27 \, {\left (\cos \left (d x + c\right )^{4} + \cos \left (d x + c\right )^{3}\right )} \sqrt {a} \log \left (\frac {a \cos \left (d x + c\right )^{3} - 7 \, a \cos \left (d x + c\right )^{2} + 4 \, \sqrt {a \cos \left (d x + c\right ) + a} \sqrt {a} {\left (\cos \left (d x + c\right ) - 2\right )} \sin \left (d x + c\right ) + 8 \, a}{\cos \left (d x + c\right )^{3} + \cos \left (d x + c\right )^{2}}\right ) + 4 \, \sqrt {a \cos \left (d x + c\right ) + a} {\left (21 \, \cos \left (d x + c\right )^{2} - 2 \, \cos \left (d x + c\right ) + 8\right )} \sin \left (d x + c\right ) + \frac {48 \, \sqrt {2} {\left (a \cos \left (d x + c\right )^{4} + a \cos \left (d x + c\right )^{3}\right )} \log \left (-\frac {\cos \left (d x + c\right )^{2} - \frac {2 \, \sqrt {2} \sqrt {a \cos \left (d x + c\right ) + a} \sin \left (d x + c\right )}{\sqrt {a}} - 2 \, \cos \left (d x + c\right ) - 3}{\cos \left (d x + c\right )^{2} + 2 \, \cos \left (d x + c\right ) + 1}\right )}{\sqrt {a}}}{96 \, {\left (a d \cos \left (d x + c\right )^{4} + a d \cos \left (d x + c\right )^{3}\right )}} \]

input
integrate(sec(d*x+c)^4/(a+a*cos(d*x+c))^(1/2),x, algorithm="fricas")
 
output
1/96*(27*(cos(d*x + c)^4 + cos(d*x + c)^3)*sqrt(a)*log((a*cos(d*x + c)^3 - 
 7*a*cos(d*x + c)^2 + 4*sqrt(a*cos(d*x + c) + a)*sqrt(a)*(cos(d*x + c) - 2 
)*sin(d*x + c) + 8*a)/(cos(d*x + c)^3 + cos(d*x + c)^2)) + 4*sqrt(a*cos(d* 
x + c) + a)*(21*cos(d*x + c)^2 - 2*cos(d*x + c) + 8)*sin(d*x + c) + 48*sqr 
t(2)*(a*cos(d*x + c)^4 + a*cos(d*x + c)^3)*log(-(cos(d*x + c)^2 - 2*sqrt(2 
)*sqrt(a*cos(d*x + c) + a)*sin(d*x + c)/sqrt(a) - 2*cos(d*x + c) - 3)/(cos 
(d*x + c)^2 + 2*cos(d*x + c) + 1))/sqrt(a))/(a*d*cos(d*x + c)^4 + a*d*cos( 
d*x + c)^3)
 
3.2.30.6 Sympy [F]

\[ \int \frac {\sec ^4(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\int \frac {\sec ^{4}{\left (c + d x \right )}}{\sqrt {a \left (\cos {\left (c + d x \right )} + 1\right )}}\, dx \]

input
integrate(sec(d*x+c)**4/(a+a*cos(d*x+c))**(1/2),x)
 
output
Integral(sec(c + d*x)**4/sqrt(a*(cos(c + d*x) + 1)), x)
 
3.2.30.7 Maxima [F(-1)]

Timed out. \[ \int \frac {\sec ^4(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\text {Timed out} \]

input
integrate(sec(d*x+c)^4/(a+a*cos(d*x+c))^(1/2),x, algorithm="maxima")
 
output
Timed out
 
3.2.30.8 Giac [F(-2)]

Exception generated. \[ \int \frac {\sec ^4(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\text {Exception raised: TypeError} \]

input
integrate(sec(d*x+c)^4/(a+a*cos(d*x+c))^(1/2),x, algorithm="giac")
 
output
Exception raised: TypeError >> an error occurred running a Giac command:IN 
PUT:sage2:=int(sage0,sageVARx):;OUTPUT:sym2poly/r2sym(const gen & e,const 
index_m & i,const vecteur & l) Error: Bad Argument Value
 
3.2.30.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\sec ^4(c+d x)}{\sqrt {a+a \cos (c+d x)}} \, dx=\int \frac {1}{{\cos \left (c+d\,x\right )}^4\,\sqrt {a+a\,\cos \left (c+d\,x\right )}} \,d x \]

input
int(1/(cos(c + d*x)^4*(a + a*cos(c + d*x))^(1/2)),x)
 
output
int(1/(cos(c + d*x)^4*(a + a*cos(c + d*x))^(1/2)), x)